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带 形状 参数 的 Bézier 曲线 的 能 量 优 化 * 


严 兰 兰 ， 殉 继 秋 ， 李 水 平 
( 东 华 理工 大 学 理学 院 , 南昌 330013) 


摘 要 : 相 较 于 经 典 的 Bézier 曲线 ， 带 形状 参数 的 BEzier 曲线 提供 了 独立 于 控制 顶点 的 形状 调整 自由 度 ， 但 同时 又 增 
加 了 设计 人 员 选 择 形 状 参数 的 工作 量 。 鉴 于 此 , 主要 讨论 形状 参数 的 选取 方案 。 首 先 证 明了 已 有 文献 中 给 出 的 Bernstein 
基 函 数 的 含 参 数 扩展 基 为 全 正 基 ， 从 而 保证 了 相应 的 带 形状 参数 的 BEzier 曲线 的 理论 价值 ; 然后 采用 能 量 最 小 化 方法 
来 确定 曲线 中 形状 参数 的 取 值 ， 推 导 了 曲线 的 拉 伸 能 量 、 弯 曲 能 量 、 扭 曲 能 量 近似 最 小 时 ， 形 状 参数 的 计算 公式 ， 为 
曲线 的 应 用 提供 了 方便 。 

关键 词 : Bézier 曲线 ; 形状 参数 ; 能 量 优 化 ; 参数 选择 

中 图 分 类 号 : TP391.72 doi: 10.3969/j.issn.1001-3695.2017.10.1018 


Energy optimization of Bézier curves with shape parameter 


Yan Lanlan, Fan Jiqiu, Li Shuiping 
(College of Science, East China University of Technology, Nanchang 330013, China) 


Abstract: Compared with the classical Bézier curves, the Bézier curves with shape parameter provide the shape adjustment 
freedom which independent of the control points. However, the introduction of shape parameter increases the workload of 
designers for choosing shape parameters. In view of this, this paper mainly discussed the selection of shape parameter. Firstly, 
it proved that the extended basis of Bernstein basis functions with parameter presented in the existing literature was totally 
positive basis. This guarantee the theory value of the corresponding Bézier curves with shape parameter. Then it used the energy 


minimization method to determine the shape parameter of the curves. It deduced the calculation formula of the shape parameter 


which made the stretch energy, strain energy and jerk energy of the curves approximate minimum. The formula was convenient 


for the application of the curves. 
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的 关键 曲线 是 否 光 顺 以 及 曲面 的 曲率 ( 主 曲率 、 高 斯 曲率 、 平 均 


0 引言 


率 等 ) 变 化 是 否 均匀 来 判断 曲面 是 否 光 顺 。 文 献 [3~6] 研 究 了 

在 计算 机 辅助 几何 设计 (CAGD) 领 域 , 曲线 曲面 造型 是 一 个 ” 能 量 最 优化 方法 在 曲线 曲面 光 顺 处 理 中 的 应 用 。 
最 基本 最 核心 的 问题 。Bézier 方法 由 于 采用 独特 的 Bernstein 多 能 量 优化 造型 方法 的 基本 思想 ， 是 以 数学 规划 和 优化 问题 
项 式 作为 基 函 数 ， 使 得 它 具 有 许多 优良 的 性 质 ， 一 经 问世 ， 就 ” 为 表达 形式 ， 以 曲线 曲面 具有 最 小 变形 能 量 为 目标 ， 运 用 各 种 
受到 工业 界 和 CAGD 学 术 界 的 广泛 重视 由， 从 而 成 为 曲线 曲 盏 约束 以 及 施加 外 荷载 的 方式 来 控制 曲线 曲面 的 形状 。 采 用 能 量 
造型 中 广 为 应 用 的 基本 方法 和 工具 之 一 。 天 化 方法 进行 曲线 设计 ， 可 以 根据 不 同 的 设计 要 求 选择 不 同 的 
Bézier 方法 是 由 控制 顶点 定义 曲线 曲面 的 方法 ， 当 选 定 上 能 量 约束 函数 ,这 种 方法 较 之 传统 的 几何 方法 ,具有 可 控 性 好 、 
线 曲 面 的 次 数 以 后 ， 曲 线 曲 面 的 形状 便 由 控制 顶点 唯一 确定 。 真实 性 强 的 特征 。 关 于 能 量 最 小 曲线 造型 方法 的 研究 文献 有 很 


区 状 ， 控 制 顶点 的 选择 至 关 重 要 。 文 种 


p= 


因此 ， 要 想 获 得 预期 的 ] 多 ， 如 文献 [7~13], 文献 [14] 则 研究 了 曲线 形状 与 曲线 能 量 之 间 
[2] 研 究 了 在 给 定 部 分 控制 顶点 的 情况 下 ， 如 何 构造 其 他 的 控制 。” 的 关系 。 

顶点 ， 使 最 终 定 义 的 Bézier 曲线 具有 最 小 的 能 量 。 优 化 某 一 能 传统 BEzier 曲线 不 存在 独立 于 控制 顶点 的 形状 调整 自由 度 ， 
量 函数 ， 即 能 量 优化 法 ， 是 构造 光 顺 曲线 的 常用 方法 。 该 方法 ”很 多 文献 如 [15~21] 针 对 这 一 问题 进行 研究 , 给 出 了 多 种 含 形状 
在 动画 设计 、 计 算 机 视觉 、 图 像 处理 、 工 业 设 计 等 领域 有 着 广 参数 ,性 质 类 似 于 Bézier 曲线 ， 且 以 Bézier 曲线 为 特例 的 新 
泛 的 应 用 。 例 如 在 船体 、 汽 车 车 身 的 设计 中 ， 通 常 根据 曲面 上 线 模 型 。 这 类 含 形状 参数 的 Bézier 曲线 模型 的 优点 ， 在 于 可 以 
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录用 稿 


在 不 改变 控制 项 点 的 前 提 下 ， 通 过 改变 形状 参数 的 取 值 来 调整 
线形 状 。 文 献 [2] 以 传统 Bézier 曲线 为 研究 对 象 ， 讨 论 如 何在 
已 知 一 部 分 控制 顶点 的 条 件 下 ， 根 据 选 定 的 能 量 优 化 目标 计算 
未 知 控 制 项 点 。 注 意 到 有 些 时 候 曲 线 控制 顶点 是 事先 给 定 的 ， 
且 当 控制 顶点 是 取 自 实物 的 精确 测量 点 时 ， 不 宜 对 其 进行 调整 
或 者 由 计算 得 到 ， 这 时 含 形状 参数 的 曲线 模型 的 优越 性 便 体 现 
出 来 了 。 在 这 种 情况 下 ， 要 想 获得 预期 的 形状 ， 可 以 根据 具体 
的 设计 要 求 选择 相应 的 能 量 函 数 ， 然 后 用 优化 方法 求 出 使 能 量 
函数 达到 最 小 值 的 形状 参数 取 值 ， 这 就 是 本 文 的 研究 内 容 。 

本 文 以 文献 [17] 中 含 形状 参数 的 曲线 模型 为 研究 对 象 ， 
择 文 献 [2] 中 给 出 的 三 种 能 量 为 衡量 标准 , 给 出 使 这 三 种 能 
最 小 值 的 形状 参数 计算 公式 


1 ， 基 函数 的 表示 
文献 [15] 给 


NY 


出 了 由 三 个 含 参数 4 的 多 项 式 函 数 构成 的 函数 


BOD)=0-A1-7) 
b=2+NDI-Dt, t 
b(t)=(1—A4+ADr 


e[0,1] (1) 


当 旬 =0 时 ， 函 数组 (1) 即 为 二 次 Bernstein 基 函 数 。 
对 n>3， 借助 递 推 公式 
b’ (0) =o (0) +tb"7 (1), t el[0,1] 

其 中 : i = 0,1,.…,n，, 并 规定 5 (7)=B" (1)=0。 文献 [17] 对 式 
(1) 所 给 函数 组 进行 了 扩展 。 当 =0 时 ， 函 数组 (2) 即 为 n 次 
Bernstein 基 函 数 。 

式 (1)(2) 共 同 给 出 了 n(n>2) 次 Bernstein 基 函 数 的 含 
扩展 。 为 了 方便 ， 在 不 至 于 引起 混淆 时 ， 下 文中 将 省 略 自 


2) 


根据 式 (3)， 
的 线性 组 合 ， 即 


有 
及 = 


GCC- prr 4 
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pr” 也 可 以 表示 成 


状 参 数 的 BE6zier 全 傣 朋 


个 有 +1 次 Bernstein 


i+l 
Cn 二 CoH 


全 02DB +gn,i)B 


[ou 
le 
洋 


G+(C 到 +Cm- EC A B"™*! 


itl 


(5) 


根据 式 (5)，n 阶 4-B 基 和 n+1 次 Bernstein 基 函 数 之 间 的 
关系 可 以 用 矩阵 表示 为 
b=BJ (0) 
其 中 : 
bp=(by br bp’) 
B (B2” B""! | 
1 0 0 0 0 
生 寻 0 0 各 
n(n-3)4 C2 +(n-3)4 
i 可 
本 C2+(2n-4-C2 1)4 
J=|0 0 可 
C34 0 
C2 
0 0 0 党 ”其 
0 0 0 0 1 
这 是 一 个 (n+2)x(n+Dl) 的 双 对 角 矩 阵 。 
2 基 范 数 的 全 正 性 
文献 [17] 和 [20] 均 未 对 4-B 基 的 全 正 性 进行 讨论 ,下 面 先 给 
出 与 全 正 基 有 关 的 概念 和 结论 ， 再 来 证 明 4-B 基 为 全 正 基 。 
定义 1 全 正和 矩阵 。 若 矩阵 五 的 所 有 子 式 都 非 负 ， 则 称 之 
为 全 正和 矩阵 。 


的 记号 ， 如 将 br (1) 简 记 为 b”。 
文献 [20] 给 出 了 含 两 个 参数 w 条 的 初始 函数 组 ， 当 
w=B 人 4 时 ， 该 函数 组 与 式 (0) 一 致 。 文 献 [20] 同 样 借助 式 (2) 


1 


对 初始 函数 组 进行 了 扩展 ， 


对 称 性 。 


根据 文献 [20]， 当 


函数 组 {Br} ， 具有 规范 性 、 


4e(-2,1] 时 
为 了 方便 ， 
根据 文献 [20]，7 


时 , 该 函数 旨 线性 无 关 , 因此 构成 
你 该 函数 组 为 怀 阶 4-B 基 。 
阶 4-B 基 可 以 统一 地 显 式 表示 成 


br =[C, + (Cs 一 C 2)4+ 


3 

ee 0 - 

其 中 : i=0,1,.….,n， n>2。 将 式 (3) 改 写 , 可 以 将 br 表示 成 一 

个 nn 次 Bernstein 基 函 数 和 一 个 n++1 次 Bernstein 基 函 数 的 线性 
组 合 ， 即 


(i+D(2i-n) AB"™!! 


(nt+Dn i+l 


1 太 =[1+ i(n—2i+1) A]B"+ 


n(n-l) 


(4) 


函数 组 {bp"} 则 为 扩展 结果 的 特例 。 


组 基 函 数 。 


定义 2 全 正 基 ,。 设 {460,44,….， 


的 基 函 数组 ， 若 对 于 任意 的 节点 系列 QS 加 < 在 <.…< 
该 基 函 数组 的 配 


均 为 全 正 秆 阵 


结论 1 全 正和 矩阵 的 乘积 仍 为 全 正和 矩阵 。 


矩阵 : 


结论 2 
命题 1 当 
素 非 负 。 


证 明 运用 数学 归纳 法 。 


经 典 的 Bernstein 


.sh, 
|- (uj(£,)) 0 jo 


T 
T 
I 


7 则 称 人 0 


0 为 全 


翘 函数 为 全 正 


当 n=2 时 ， 

1 0 0 

L424 244 
J= 3 3 

0 0 1 


4,} 为 定义 在 闭 区 间 [a,b] 上 


tb, 


Xe(-2,1 吕 时， 对 所 有 n>2， 和 矩阵 J 的 所 有 元 
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在 条 件 和 4e(-2,1] 下 , 该 矩阵 所 有 元 素 非 负 。 假设 当 n = 
时 ,命题 成 立 , 即 对 n= ，i=0,1,…k， 式 (5) 右 端的 系数 总 
是 非 负 ， 则 当 n=k+1 时 ， 结 合式 (2)(5) 可 以 推出 : 


g (KC +g (Kil) Ch Be 


CH 


bt+! = f KDCEn+f Ki-DC B+ 


Ch 


在 假设 之 下 ， 对 所 有 i=0,1,...,k+1 ， 上 式 右 端的 系数 非 
负 。 证 毕 。 
命题 2 当 4e(-2,1 时 ， 对 所 有 n>2，J 为 全 正和 矩阵 。 
证 明 根据 命题 1，J 的 所 有 元 素 非 负 , 由 于 J 为 双 对 角 和 天 
阵 ， 易 知 其 所 有 子 式 非 负 ， 所 以 ./ 为 全 正和 矩阵 。 证 毕 。 
命题 3 当 4e(-2,1] 时 , 对 所 有 n>2, 7n 阶 4-B 基 为 全 正 


区 


， 


证 明 任 给 节点 系列 0<th <t <...<t <1， 分 别 用 M,、 
Ms 表示 nn 阶 4-B 基 、n+1 次 Bernstein 基 函 数 的 配置 矩阵 ， 


则 根据 式 (6) 有 
M,=M,] 

于 Ms 和 J 均 为 全 正 箱 阵 ， 所 以 

正 矩 阵 ， 故 4-8 基 为 全 正 基 。 


3 ”形状 参数 的 选择 


给 定 控制 顶点 P eR” (d =2,3;i=0,],...， 
条 nL 阶 4-B 曲线 如 下 : 


结论 1 知 用, 亦 为 全 


nn) ,可 以 定义 一 


p(n)= >.bP (7) 
让 
其 中 : te[0,1]; 4e(-2,1]。 
4-B 基 的 性 质 可 知 : 4-B 曲线 具有 类 似 于 Bezier 曲线 的 


仿 射 不 变性 、 变 差 缩 减 性 。 除 
线 还 具有 形状 可 


性、 对 称 性 、 几 何不 变性 、 
> 外， 由 于 4-B 基 中 含 参数 4 ， 所 以 4- 


选择 文献 [2] 中 使 用 的 三 种 能 量 函 数 : 


已 = pw OF dt, k=1,2,3 


当 大 = 工时 玉 , 为 拉 伸 能 量 的 近似 ， 其 反映 的 是 曲线 长 度 ; 
当 k=2 时 ,为 弯曲 能 量 的 近似 ,其 反映 的 是 曲线 的 曲率 ， 当 
k=3 时 为 扭曲 能 量 的 近似 ,其 反映 的 是 曲线 曲率 的 变化 量 ， 
下 面 推导 使 (k =1,2,3) 取 最 小 值 的 参数 4 的 求解 公式 。 


将 式 (4) 代 入 式 (7) 并 整理 ， 可 得 


PO=> 


对 上 式 求 k(k==1,2,3) 阶 导 
扎 三 i(n—2i+l 
P = WA {+ 和 二]P)B” 二 


(n+1)! k (2i—n—2) B” —k+l 
(n—k+1)! > A | (n+l)n 4P, ,|B; 


n+ 
i(n—2i+l1) n i(2i—n—2) nt+l 
+ APB, + PUB 


数 ， 可 得 


(8) 


记 
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“J 
:Tt 曲线 的 能 量 优化 


Q, =i(n—2i+ DP 


W, =i(2i—n -2)P, 
f= A SAPpB” k (9) 
i=0 


n—k+l 
_ (n-2)! k n— 5 (n-D)! k n—k+l 
t= 济 光 4 Q.B, rp k+1)! > A WB 
i=0 


则 式 (8) 可 整理 成 
pWD=f+g4 
进而 p "Dp OD)=f:f+24f :gt+hg:g 
则 曲线 能 量 为 关于 形状 参数 4 的 函数 ， 即 


1 . 2 1 
E(D=| ffat22| f -gait gs 


令 当中 = 站 sw+24| gs =0 
"fgdt 

4=_h (10) 
J gga 


其 中 : f 和 8 由 式 (9) 定 义 。 


为 了 方便 使 用 ， 下 面 分 别 给 出 当 k=1,2,3 时 ， 取 
n= 2 3,4,5 所 得 式 (10) 的 体 表 达 式 。 
先 统一 记 
D,=(P, P P,) 


其 中 : n=2,3,4,5。 
4 天 =1 时 ， 取 n=2， 式 (10) 即 为 


I 


1 
1 -4 2 
XI=|0 4 -4 
0 0 1 
人 2 -5 1 
Yi=I0 5 -5 
0 0 2 


当 庆 =1 时 ， 取 n=3， 式 (10) 即 为 


DIX!D,; 
Le 
2 -5 -1 2 
| 10 3 0 -1 
XX! = 
3 |0 0 3 -5 
0 0 0 2 
3 -6 1 -1 
1 10 3 -1 1 
Y= 
3 |0 0 3 -6 
0 0 0 3 


当 庆 =1 时 ， 取 n=4， 式 (10) 即 为 


Tyl 

1 9 Di XaDs 
Tyl 

2 DiYiD, 


~ \V 
CNAINa 
1 LA、 


\ OO/ 


i I 
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录用 稿 优化 


当 k=2 时 ， 取 n=4， 式 (10) 即 为 


2-_2 DI XD, 


DIYiD, 


15 -30 -10 4 6 4 = 7 
2 DIY?D, 
0 16 2 -8 4 
XI=-|0 0 8 2 -10 其 中 : 
0 0 0 16 -30 28 -81 20 9 -4 
0 0 0 0 15 
过 二 0 56 -32 -8 9 
20 -35 -3 5 -7 X2=-|0 0 12 -32 20 
0 17 -3 -1 5 0 0 0 56 -8l 
Yi=|0 0 6 -3 -3 0 0 0 0 28 
0 0 0 17 335 52 -117 21 -2 -6 
0 0 0 0 20 0 68 -28 11 -2 
当 k=1 时 ， 取 n=5 ， 式 (10) 即 为 oly 
a 0 0 0 68 -117 
i $5 LsAsUs 
1 = -2 Dip 0 0 0 0 52 
28 49 19 1 7 4 当 k =2 时 ， 取 n=5， 式 (10) 即 为 
0 25 5 -1l -2 7 1 = 5 OD 
2 Diy2 
wi 2 
5 |0 0 0 10 5 -19 其 中 : 
0 0 0 0 5 -4 125 -360 80 36 4 -10 
i TT 0 259 -140 -40 18 4 
175 -280 -70 40 -5 -35 0 0 36 -8 -40 36 
0 145 -10 -32 37 -5 5 | 0 0 0 36 -140 80 
i 0 0 40 -8 -32 40 0 0 0 0 259 -360 
0 0 0 40 -10 -70 0 0 0 0 0 125 
0 0 0 0 145 -280 125 -300 50 10 -10 0 
0 0 0 0 0 175 0 191 -82 -2 12 -10 
Ss 
5 |0 0 0 16 -82 50 
0 0 0 0 191 -300 
l 0 0 0 0 0 125 
0 4 当 久 =3 时 ， 取 n=2， 由 式 (10) 得 出 和 =0， 此 时 的 2 阶 
和 0 0 1 4.-B 曲线 即 为 2 次 Bézier 曲线 。 
9 1 -1 1 当 k =3 时 ， 取 n= 二 3 ， 式 (10) 即 为 
了 二 0 1 一 1 DIxX3D. 
0 0 1 14= 3 an 
当 k=2 时 ， 取 n=3， 式 (10) 即 为 其 中 : 
1 _5 DX3D, 1 -4 4 一 2 
”PD ; |0 3 -6 4 
用 3 = 
] -3 1 0 0 0 3 -4 
i 0 2 -2 1 0 0 0 1 
3 |0 0 2 -3 元 4 光 
0 0 0 1 yl0 7 2 2 
大 中 rd | Wi 
y2 10 2 -1 0 0 0 7 
”10 0 2 -4 当天 =3 时， 取 寻 = 4 ， 式 (10) 即 为 
0 0 0 2 


:201805.00014v1 


chinaXiv 


录用 稿 
其 中 : 
7 -29 20 -3 -2 
0 26 -32 12 -3 
X=I0 0 12 -32 20 
0 0 0 26 -29 
0 0 0 0 7 
16 -44 12 -4 4 
0 31 -16 2 -4 
Y=I0 0 4 -16 12 
0 0 0 31 -44 
0 0 0 0 16 
当 k=3 时 ， 取 n=5， 式 (10) 即 为 
DIX3D; 
4=— 
其 中 : 
4 -17 11 -1 -1 0 
0 17 -21 3 2 -1 
，|10 0 6 -4 3 =- 
式 三 
0 0 0 6 -21 11 
0 0 0 0 17 -17 
0 0 0 0 0 4 
55 -176 66 -4 -1l 5 
0 145 -ll4 12 -ll1 -1 
a 0 0 24 -8 12 -4 
5 |0 0 0 24 -114 66 
0 0 0 0 145 -176 
0 0 0 0 0 55 
4 数值 实验 

首先 以 一 段 由 4 个 控制 顶点 定义 的 3 阶 4-B 开 曲线 和 3 次 
Bezier 开 曲 线 为 例 ， 分 析 由 相同 控制 顶点 和 不 同 能 量 目标 定义 
的 4-B 曲线 与 Bézier 曲线 的 差异 。 

图 1 所 示 红 、 绿 、 蓝 色 曲 线 分 别 为 当 k =1、k=2、k=3 
时 的 4-B 曲线 ( 见 电 子 版 )， 公式 计算 出 的 参数 分 别 为 
4=- 下 、4=- 站 、4 = 站; 黑色 曲线 为 Bézier 曲线 。 由 图 
1 可 知 ， 以 不 同 能 量 函数 为 优化 目标 得 到 的 4-B 曲线 具有 不 同 
的 对 控制 多 边 形 的 逼近 程度 ;大 越 大 , 曲线 越 接近 控制 多 边 形 ; 
当 k =3 时 的 4-B 曲线 与 Bézier 曲线 非常 接近 。 

图 1 3 阶 4-B 曲线 段 与 3 次 Bézier 曲线 段 


图 2 中 子 图 2(a)~(@ 分 别 为 图 1 中 红 、 绿 、 蓝 色 曲 线 的 曲率 
图 ; (qd) 为 图 1 中 黑色 曲线 的 曲率 图 。 由 图 2 可 知 ， 以 拉 伸 能 量 
近似 最 小 为 目标 构造 的 曲线 具有 较 大 的 曲率 ， 曲 率 的 最 大 值 在 
端点 处 取得 ; 以 弯曲 能 量 近 似 最 小 为 目标 构造 的 曲线 具有 相对 
较 小 的 曲率 ， 且 曲线 的 曲率 变化 范围 最 小 ， 以 扭曲 能 量 近似 最 
小 为 目标 构造 的 曲线 在 端点 处 出 现 零 曲 率 ， 其 具有 与 Bézier 
线 相 似 的 曲率 图 。 
8 1 
6 0.8 
4 06 
2 0.4 
4 0 0.5 1 


(b)k=2 


(co) k=3 


图 2 


怠 上 县 - 


注 : 虽然 扭 映 


图 1 中 各 


(d)Bézier 
线段 的 曲率 图 


有 星 


的 是 


率 的 变化 量 , K=3 时 的 


E, 为 扭曲 能 量 的 近似 ， 但 
线 的 曲率 变化 量 并 不 比 大 
导致 这 种 结果 的 主要 原因 是 El 


从 图 2 
=2 时 的 4-B 曲线 的 曲率 变化 量 小 ， 


可 以 看 出 ，k=3 时 的 4-B 曲 


为 能 量 的 近似 值 而 非 精 确 值 。 


为 了 更 加 直观 地 显示 采 / 


能 量 法 确定 


状 参数 的 造型 效果 ， 


通过 选择 合适 的 控 人 


吊顶 点 , 下面 给 出 整体 G 连续 的 2~5 阶 4-B 


开 曲 线 、 闭 曲线 ， 


以 及 对 应 的 


相同 控制 顶点 定义 的 2~5 次 


Bezier 曲线 ， 其 中 4-B 曲线 的 
确定 ， 取 4=0 即 得 Bézier 出 


图 3 所 示 为 由 相同 开 控制 
的 2 阶 4-B 曲线 和 2 次 Bézier 
数 依次 为 一 1.3036 、 
—1.73 、--0.5331 、-0.6957 、 
一 1.3950 、 一 1.1872 、 
-0.5448 、-0.5634 、-0.2624 
-0.1420 、-0.1928 、--0.0725 
—0.4521 、--0.1737 ; 

图 4 所 示 为 由 相同 开 控制 
的 3 阶 4-B 曲线 和 3 次 Bézier 
数 依次 为 -4.3380 、0.3796、 
一 3.6123 、 一 3.7917 ; (b) 中 各 
1.4751、 —0.6825 、 0.3269 


一 1.32435 、 
—0.3166 、 
—0.65 ; 


(o) 中 的 4-B 


—2.9355 、 
条 曲线 段 的 参数 依次 为 -1.8555 、 
\、-0.9821 、-1.3004 、 


的 各 条 曲线 段 的 参数 4 根据 式 (10) 


线 。 
图 3~10 中 4-B 曲线 的 起 始 段 均 用 第 


头 标 出 。 
多 边 形 按 不 同 能 量 优 化 目标 定义 
线 。 图 3(a) 中 各 条 曲线 段 的 参 
—0.8605 、-0.14206 、-0.3333 、 
—1.4769 、—1.1662 、 
(b) 中 各 条 曲线 段 的 参数 依次 为 
、—0.0302 、-0.0823 、 一 1.1234 、 
、—0.7250 、-0.4372 、-0.6312 、 
线 即 为 Bézier 曲线 。 
多 边 形 按 不 同 能 量 优 化 目标 定义 
线 。 图 4(a) 中 各 条 曲线 段 的 参 
—1.6767 、 -3.1111、 


一 1.4773 ; 
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(c) 中 各 条 曲线 段 的 参数 依次 为 -0.0928 、1.4437 、0.3368 、 0.2971、-2.4821 ，(c) 中 各 条 曲线 段 的 参数 依次 为 -0.3635 、 
0.2712 、0.1544 、0.0755 、-9.6892x10 5 ; (c) 为 3 次 Bézier 一 1].1717 、0.7754 、2.4845 、-0.0324 ; (d) 为 5 次 Bézier 曲 


GaK=1 (k=2 (Ok =3 (Btzien 
图 3 2 阶 4-B 1 线 与 2 次 Bézier 线 


(co) k=3 (d)Bézier 
图 5 4 阶 4-B 开 曲 线 与 4 次 Bézier 开 曲 线 


(c)k=3 (d)Bézier 
图 4 3 阶 4-B 开 曲 线 与 3 次 Bezier 开 曲 线 


5 所 示 为 由 相同 开 控 制 多 边 形 按 不 同 能 量 优 化 目标 定义 
的 4 阶 4-B 曲线 和 4 次 Bézier 曲线 。 图 5(a) 中 各 条 曲线 段 的 参 
数 依次 为 -6.1598 、4.7403 、-6.2274 、1.8876 、--0.9417 、 

一 0.4019 、-4.0828 ; (b) 中 各 条 曲线 段 的 参数 依次 为 -2.1410 、 
2.9145 、-1.8178 、1.4827 、0.85$1、0.9537 、--0.8249 ; 

子 (c) 中 各 条 曲线 段 的 参数 依次 为 -0.1199 、-0.7964 、0.3131、 
-0.8238 、-0.3056 、0.6755 、0.3145 ; 子 (d) 为 4 次 Bézier 


线 。 (ck=3 (d)Bézier 
图 6 5 阶 4-B 开 曲 线 与 5 次 Bézier 开 曲 线 
图 6 所 示 为 由 相同 开 控 制 多 边 形 按 不 同 能 量 优 化 目标 定义 
的 5 阶 4-B 曲线 和 5 次 Bézier 曲线 。 图 6(a) 中 各 条 曲线 段 的 参 图 7 所 示 为 由 相同 闭 控制 多 边 形 按 不 同 能 量 优化 目标 定义 
数 依次 为 -6.1575 、2.4708、1.2320 、-5.1501、-8.2745 ;的 2 阶 4.B 曲线 和 2 次 Bézier 曲线 。 图 7(a) 中 各 条 曲线 段 的 参 


(b) 中 各 条 曲线 段 的 参数 依次 为 -1.9719 、--1.2823 、0.5281、 数 依次 为 -0.8695 、-0.8902 、 一 1.2725 、 一 1.1717 、-0.9872 、 
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—0.5223 、—1.1072 、—1.1072 、—0.5223 、 -0.9872 、-1.1717 、 


一 1.2725 、--0.8902 、-0.8695 ; (b) 中 各 条 曲线 段 的 参数 依次 
为 -0.2666 、-0.2781 、-0.5183 、-0.4411 、 -0.3263 、 


—0.1320 、 -0.3975 、-0.3973 、-0.1320 、-0.3263 、-0.4411 、 
一 0.5183 、-0.2781 、--0.2666 ; (oc) 中 的 4-B 曲线 即 为 Bézier 


(c) k=3 (Bézier) 
图 7 2 阶 4-B 闭 曲 线 与 2 次 Bézier 闭 曲 线 


图 8 所 示 为 由 相同 闭 控制 多 边 形 按 不 同 能 量 优 化 目标 定义 
的 3 阶 4-B 曲线 和 3 次 Bézier 曲线 。 图 8(a) 中 各 条 曲线 段 的 参 
数 依次 为 一 2.0095 、--0.9397 、0.6806 、-2.8283 、0.6806 、 
一 0.9397 、 一 2.0095 ; (b) 中 各 条 曲线 段 的 参数 依次 为 -0.6110、 
—0.8906 、0.8144 、-0.2083 、0.8144 、 0.8906 、-0.6110;; 
(0) 中 各 条 曲线 段 的 参数 依次 为 0.1133 、-0.8485 、0.8392 、 


0 


-三 ”0.6311、0.8392 、-0.8485 、0.1133 ; (d) 为 3 次 Btzier 曲线 。 


(d)Bézier 


图 8 3 阶 4-B 闭 曲 线 与 3 次 Bézier 闭 曲线 
图 9 所 示 为 由 相同 闭 控 制 多 边 形 按 不 同 能 量 优化 目标 定义 
的 4 阶 4-B 曲线 和 4 次 Bézier 曲线 。 图 9(a) 中 各 条 曲线 段 的 参 
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数 依次 为 -4.6098 、 一 3.6230 、-0.6415 、 一 3.6230 、 一 4.6098 ; 
(b) 中 各 曲线 段 的 参数 依次 为 -1.6801 、-1.3301 、--0.7705 、 
-1.3301、-1.6801 ; (co 中 各 条 曲线 段 的 参数 依次 为 -0.1815 、 


-0.3343 、-0.7165 、--0.35343 、 一 0.1815 ; (d) 为 4 次 Bezier 


(d)Bézier 
线 与 4 次 Bézier 闭 曲线 


图 10 所 示 为 由 相同 闭 控制 多 边 形 按 不 同 能 量 优 化 目标 定 
义 的 5 阶 4-B 曲线 和 5 次 Bézier 曲线 。 图 10(a) 中 各 曲线 段 的 
参数 依次 为 -4.8531 、 一 7.1176 、 一 4.8531 ; (b) 中 各 条 曲线 段 
的 参数 依次 为 -2.4292 、-3.2692 、 一 2.4292 ; (c) 中 各 条 曲线 
段 的 参数 依次 为 -1.4611 、-0.6604 、 一 1.4611; (d 为 5 次 
Bezier 曲线 。 


(k=1 (bk=2 


(ck=3 


(d)Bézier 


图 10 5 阶 4-B 闭 


线 与 5 次 Bézier 闭 曲 线 


从 图 3~10 可 以 看 出 ， 取 拉 伸 能 量 近 似 最 小 为 目标 得 到 的 
线 “ 棱 角 分 明 ”， 各 条 曲线 段 呈 拉 直 状 ， 在 分 段 连 接点 处 容易 
形成 “ 尖 角 ”*”， 曲 线 的 阶 数 7 越 低 ， 这 种 特征 越 明 显 ， 取 弯曲 和 


CC 


量 近似 最 小 为 目标 得 到 的 曲线 总 是 具有 良好 的 视觉 效果 ， 取 捍 
能 量 近似 最 小 为 目标 得 到 的 曲线 大 部 分 时 候 与 相应 的 Bézier 
线 具 有 相似 的 外 形 ， 从 计算 出 的 参数 数据 来 看 ， 它 们 在 某 些 
段 上 与 Bezier 曲线 4=0 的 取 值 较 接近 ， 但 在 有 些 段 处 却 存在 
较 大 差异 。 
在 前 面 提 到 , 当 4e(-2,1 时 ，4-B 基 非 负 且 全 正 , 基 函 数 
的 非 负 性 可 以 保证 曲线 的 凸 包 性 。 而 由 能 量 优化 方法 推导 出 的 
参数 计算 公式 与 控制 顶点 的 坐标 密切 相关 ， 无 法 保证 所 有 情况 
下 得 到 的 参数 4 均 落 在 区 间 (-2 直 内 ， 但 从 数值 实例 可 以 看 出 ， 
这 并 未 影响 曲线 的 凸 包 性 。 究 其 原因 ， 基 函数 非 负 是 
凸 包 性 的 充分 而 非 必要 条 件 。 
注 : 在 实际 应 用 中 ， 可 以 让 不 同 的 曲线 段 按 不 同 的 能 量 优 
化 目标 来 确定 形状 参数 ， 也 可 以 将 不 同 的 能 量 优化 目标 加 权 组 
合 来 确定 形状 参数 ， 设 计 者 可 以 根据 具体 的 设计 要 求 选择 最 合 
适 的 形状 参数 确定 方案 。 
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本 文 以 含 形状 参数 的 Bézier 曲线 模型 为 研究 对 象 ， 首 先 从 
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理论 上 证 明了 该 模型 的 应 用 价值 ;然后 给 出 了 分 别 以 拉 伸 能 量 、 
弯曲 能 量 、 扭 曲 能 量 近似 最 小 为 优化 目标 得 到 的 形状 参数 计算 
公式 ， 通 过 曲线 图 和 曲率 图 直观 比较 了 不 同 能 量 优 化 目标 的 差 
异 ， 并 给 出 了 日 常生 活 中 一 些 常 见 图 案 的 造型 模拟 。 

近 十 五 年 左右 , 带 形状 参数 的 Bézier 曲线 成 为 CAGD 中 的 
开 究 热点 之 一 。 除 了 本 文 引用 的 文献 [15~21] 以 外 ,还 有 一 大 批 
文献 给 出 了 性 质 类 似 于 Bézier 曲线 且 含 形状 参数 的 新 曲线 。 与 
现 有 文献 相 比 ， 本 文 优点 主要 体现 在 两 个 方面 : 一 是 考察 了 基 
函数 的 全 正 性 ， 基 函数 的 全 正 性 决定 着 相应 曲线 的 变 差 缩减 忻 
和 保 晤 性， 因此 决定 着 曲线 的 应 用 价值 。 从 这 个 意义 上 讲 ， 考 
察 基 函数 是 否 具备 全 正 性 是 衡量 新 曲线 是 否 有 其 存在 价值 的 标 
准 之 一 。 二 是 给 出 了 基于 能 量 优 化 目标 的 形状 参数 取 值 的 计算 
公式 。 在 曲线 中 引入 形状 参数 可 以 赋予 曲线 独立 于 控制 顶点 的 
形状 调整 能 力 ， 但 同时 也 增加 了 设计 人 员 选 择 形状 参数 的 工作 
量 。 提 供 以 常见 设计 要 求 为 目标 所 得 到 的 形状 参数 计算 公式 ， 

有 助 于 设计 人 员 根 据 具体 设计 目标 快速 确定 合适 的 形状 参数 ， 

从 而 为 曲线 的 应 用 提供 了 便利 。 
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